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ABSTRACT 

Aims. The high energy spectrum of 3C 273 is usually understood in terms of inverse-Compton emission in a relativistic leptonic jet. 
This model predicts variability patterns and delays that could be tested with simultaneous observations from the radio to the GeV 
range. 

Methods. The instruments IBIS, SPI, JEM-X on board INTEGRAL, PGA on board RXTE, and FAT on board Fermi have enough 
sensitivity to follow the spectral variability of 3C 273 from the keV to the GeV. We looked for correlations between the different energy 
bands, including radio data at 37 GHz collected at the Metsahovi Radio Observatory and built quasi-simultaneous multiwavelength 
spectra in the high energy domain when the source is flaring either in the X-rays or in the y rays. 

Results. Both temporal and spectral analysis suggest a two-component model to explain the complete high energy spectrum. X-ray 
emission is likely dominated by a Seyfert-like component while the y-ray emission is dominated by a blazar-like component produced 
by the relativistic jet. The variability of the blazar-like component is discussed, comparing the spectral parameters in the two different 
spectral states. Changes of the electron Lorentz factor are found to be the most likely source of the observed variability. 
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1. Introduction 

Active galactic nuclei (AGN) are very energetic extragalac- 
tic sources powered by accretion on supermassive black holes. 
Blazars are a class of AGN characterized by a relativistic jet 
pointing towards the observer. Because of the beaming, the 
emission originating in the jet dominates the observed spectrum. 

The quasar 3C 273 is the brightest and hence one of the best 
monitored AGN. Located at a distance of z = 0.158 ( [Strauss ^ 
ar]|1992| l, it features interesting blazar-like (a jet with supeiTu- 
minal motion and high variability) and Seyfert-like (the strong 
blue bump, the soft X-ray excess, and variable emission lines) 
characteristics. 

Several observational campaigns have been carried out on 
3C 273, starting from the broad multiwavelength campaign by 
Lichti et ah] (|1995|l, which made use of data spanning from ra 


dio to y rays, to more recent observations in X-rays. Cour- 
voisier et al.| ( |2003] l described simultaneous observation of IN¬ 
TEGRAL, XMM-Newton, and RXTE performed in January 


2003. jChernyakova et ah (2007 1 used simultaneous INTEGRAL 
and XMM-Newton data in the period 2003-2005 and compared 
their results to historical data to investigate the secular evolu¬ 
tion of 3C 273. Liu & Zhang| ( 201 Tj ) used XMM-Newton data to 
investigate the nature of the soft X-ray excess. 

Many works have been published to describe in detail the 
spectral and temporal characteristics of 3C 273, which has been 
widely used as a test case to validate theoretical models and 
predictions. Paltani et al. (19981 found two distinct variable 
components in the optical-ultraviolet band. Soldi et al. P008|l 


studied the temporal variability in a wide energy range. 


Grandij 


|& Palumib^ ( |2004| ) were able to decouple the non-thermal (jet) 
and thermal (accretion flow) emission by studying the long¬ 
term spectral variability of BeppoSAX data. [Marcowith et al.j 
( 1995| l tested their model of relativistic electron-positron beam 
for blazar emission. [Marscher & Gear] ( |1985] l and [Tlirler et al.j 
( |2000| l modelled the radio emission in the context of synchrotron 
emission from relativistic electrons produced by shocks in the 
jet. Artyukh et al. ( [2012 1 investigated the physics and the con¬ 
dition of the core of 3C 273 on parsec scale. A detailed review 
has been written by |Courvoisier| ( [T998] l and a database collecting 
public data of the source over several decades is available online 
( jTurler et al.|T999l[Soldi et al.|20081 l. 

The aim of this work is to use contemporaneous X-ray to y- 
ray observations obtained with INTEGRAL and Eermi for the 
first time since the CGRO campaign ( [Lichti et al.|1995[ l to study 
the origin of the non-thermal spectral components. In Sect. 
we introduce the data used for this work and explain how the 
datasets have been buildt. In Sect, [^we present the results of 
the cross-correlation analysis between lightcurves obtained with 
Eermi-LAT, RXTE-PCA, and at 37 GHz and compare our re¬ 
sults with those of previous studies. We find that y-ray and radio 
emission are correlated and that the X-ray component has a dif¬ 
ferent origin. 

In Sect. I^we describe the high energy spectrum of 3C 273 
and its variability. Data from RXTE-PCA (0.2 - 75 keV), JEM- 
X, ISGRI and SPI on board INTEGRAL (3 keV - 1 MeV), and 
Eermi-LAT (100 MeV - 100 GeV) are used. We selected two 
different spectral states based on the behaviour of the PCA and 
EAT lightcurves, and fit several models to their spectra to investi¬ 
gate the differences between the two states. We find that the best 
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Fig. 1. Top panel: RXTE-PCA lightcurve in 0.2 - 75 keV. Middle panel: Fermi-LAT lightcurve in 0.1 - 100 GeV. Area marked in light grey 
colour shows the selected flaring time epochs in both lightcurves. Bottom panel: radio lightcurve at 37 GHz. Time units are days on bottom axis 
(Modihed Julian Date) and decimal years on top axis. 


fit model is a two-component model, which can be interpreted 
as a Seyfert-like component dominating the X-ray emission, and 
blazar-like component dominating at y rays. In Sect, [^we dis¬ 
cuss the results and give a physical interpretation of the data, 
focusing on the blazar-like component, in the context of rela¬ 
tivistic jet emission and synchrotron self Compton (SSC) model, 
and finally we summarize the results in sect. 


2. The data sample 

For this work we used data from RXTE-PCA, INTEGRAL and 
Eermi-LAT, which provide a broad spectral coverage from a few 
keV up to 100 GeV and radio data collected at the Metsahovi 
Radio ObservatorjQ 

Eermi-LAT (0.1 - 100 GeV, [Atwood et ar]|2009) l data are 
used to study temporal and spectral variability at high energies. 
Thanks to its large field of view, Eermi-LAT has observed 3C 
273 almost continuously. The SED is extracted with the Max¬ 
imum Likelihood fitting technique commonly used to analyse 
Eermi data and the light curve is built using the LAT Aperture 
Photometric method assuming the average spectral index of 2.45 
as listed in the Eermi-LAT 2nd Catalog ( |Nolan et al.|2012l l and 
a binning time of 1 day. The data spans from August 2008 to 


* http://metsahovi.aalto.fi/en/ 
Article number, page 2 of 1 11 1 


July 2012 and they are extracted using the Eermi Science Tools 
version v9r27pl and the instrument response function P7_V6. 

INTEGRAL dWinkler et al.|20031l JEM-X (6 - 20 keV), IS- 
GRI (17 - 210 keV) and SPI (50 - 770 keV) data are used for the 
spectral analysis. The JEM-X and ISGRI data are analysed with 
OS A lOp] while SPI sp ectra are built with the SPI_TOOLSLIB 
tool ( |Kn5dlseder|2004| l and the spectra extracted by a model fit¬ 
ting method. They span from January 2003 to January 2012. 


RXTE-PCA (0.2 - 75 keV,||Jahoda et al. 

19961 data are col- 

lected with the HEAVENS web interfac^ ( 

Walter et al.||2010 


and span from March 2005 to December 2011. During this pe¬ 
riod 3C 273 has been observed by RXTE-PCA every few days 
for ~ 2 ks with only seven long gaps (see Pig. [^upper panel). 


The radio data at 37 GHz are obtained at the Metsahovi Ra¬ 
dio Observatory and span from September 2006 to June 2012. 
They are used for timing analysis to probe the correlation with 
the high energy emission. The flux density scale is set by ob¬ 
servations of DR 21. Sources NGC 7027, 3C 274 and 3C 84 
are used as secondary calibrators. A detailed description of the 
data reduction and analysis is given in |Terasranta et al!] ( 1998| l. 
The error estimate in the flux density includes the contribution 
from the measurement rms and the uncertainty of the absolute 
calibration. 


^ http://www.isdc.unige.ch/integral/analysis 
^ http://www.isdc.unige.ch/heavens/ 
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Fig. 2. Cross-correlation between Fermi-LAT, RXTE-PCA and radio 
lightcurves. A positive lag for an "X-Y" correlation means a delay of Y 
with respect to X. 


Soldi et aLl 


3. Temporal Analysis 

3.1. Cross-Correlation 

In the multiwavelength variability study of 3C 273 by 
( ]2008[ l, correlations between the lightcurves obtained at various 
frequencies were studied. In particular hard X-rays (> 20 keV) 
and radio-mm were found to be correlated, while no correlation 
was found with lower energy X-rays (< 20 keV) or radio-mm 
lightcurves. A difference between the hard (> 20 keV) and soft 
X-rays was also found in the fractional variability spectrum. 
This suggests a different origin for the X-rays above and below 
~ 20 keV. 


In order to study the origin of the high energy emission we 
used lightcurves obtained with Fermi-LAT, RXTE-PCA and in 
the radio with the Metsahovi Radio Telescope. The lightcurves 
are shown in Fig. 

The LAT and PCA lightcurves show a different be¬ 
haviour. The LAT data show a stable quiescent flux at 
~ 0.3 ■ 10“^ photons s ' cm“^ with small flares randomly 
distributed and three very strong flares in September 2009 (MJD 
from 55071 to 55097). These strong flares are described in 
detail by Abdo et aL] ( |2010[ ) and Rani et al. ( |2013| l. Rani et 
[ah] p013| l found a variability time scale of the order of ~ 0.1 
day for the shortest flare and of ~ 0.5 - 1 day for the others. 
The PCA lightcurve shows several consecutive flares with 
alternating high and low states, the rate varying from ~ 5 to 
~ 30 photons s ' cm“^. 


To study correlations, we used the discrete correlation func¬ 
tion (DCF, Edelson & Krolik|1988 1 with a time bin of 10 days. 
In a cross-correlation between the "X" and "Y" lightcurves, la¬ 
belled "X-Y" in Fig. a positive time lag means that the 
Y lightcurve is delayed with respect to X. Cross-correlation 
curves are shown in Fig. The solid line represents the cross¬ 



Fig. 3. Profile of the best radio response, the y axis is in arbitrary 
scale. The maximum is located at 160 days after the beginning of the 
curve. The FWHM of the curve is 197 days. See text for details. 



Fig. 4. Radio lightcurve (hlack) compared with the convolution ob¬ 
tained from the LAT lightcurve (red). Fermi-LAT lightcurve is shown in 
the bottom of the picture for comparison (not in scale). Vertical dotted 
line is MJD = 54700 days. See text for details. 


correlation between the radio and the LAT data. The correlation 
is peaked at —1-120 days with a coefficient of ~ 0.63 suggesting 
that y rays are leading the radio emission. Dashed and dotted 
lines represent the absence of significant correlation of the PCA 
lightcurve with the LAT or radio lightcurves. 

This result is in agreement with the observation of [Soldi ^ 
|ar] ( |2008| , suggesting that the hard X-ray photons come from the 
jet, whereas Softer X-rays instead are dominated by a different 
component. 

3.2. Using the 1C response 

The connection between the radio and the y-ray emission in 
radio-loud AGN has been widely investigated proving that radio 
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and y-ray luminosities are correlated (see for example Padovani 

et al. ( 1993|l; Valtaoja et al. 

1999|l; Jorstad et al. (j2001|l; Lah- 

teenmaki & Valtaoja|(|2003 1; 

Marscher et al. P008J; Hovatta et 

al.|(|2009|l;|Ackermann et al.|( 

201111). 


In order to better study the link between the GeV and the 
radio emission, we try to evaluate whether the observed radio 
variability could arise from a simple convolution of the Fermi 
lightcurve. The response function would convolve each spike 
of the LAT lightcurve to a broader and delayed radio flare. The 
sum of all these broadened flares, which may be superimposed, 
should match the observed radio lightcurve. The response would 
then naturally introduce the delay seen with the cross-correlation 
of Fig. caused by the cooling time of the electrons emitting 
synchrotron radiation and the expansion of the emitting region, 
becoming optically thin at lower energies. 

As the convolution function is not known a priori, we rely 
on the shape of a typical outburst in 3C 273 at radio frequencies. 
Tiirler et al. (1999 2000| developed a method - in the frame 
of a shock-in-jet scenario - to extract from the long-term 
multi-wavelength lightcurves the temporal evolution of a typical 
synchrotron outburst at any frequency. We use here the model 
lightcurve of an average outburst of 3C 273 at 37 GHz as derived 
from the observations obtained between 1980 and 2000 and with 
the approach and parameters described by ITiirler & Lindfords] 
( [2007| l. This model lightcurve rises relatively steeply to reach a 
maximum in less than a year before decaying gradually. 


With this response function we could not reproduce the ra¬ 
dio data because the peak of the response profile is too late and 
broad when compared to the observations. The observed radio 
lightcurve can only be reproduced by compressing the response 
by a factor of 2 in time scale and shifting it by -30 d. This com¬ 
pressed and shifted response, shown in Fig. |3 allow us to match 
the width and the position of the observed radio flares. 

The necessity to shrink the response can be explained 
considering that the convolution smears the complex y-ray 
variability, with the result that a sequence of short y-ray flares 
separated by short time lag will be observed in the radio as a 
single broad flare. As explained above, the response curve was 
constructed on long-term multi-wavelength lightcurves, which 
do not include y-ray data. It is hence possible that what was 
seen in the construction of the response as a single flare, i.e. a 
single shock wave in the modelling of |Turler et al.| ( T999l|2000| l, 
might actually be the superposition of several smaller shocks. 


It has to be pointed out that the available LAT data start later 
than the radio data. So we cannot reproduce correctly the early 
radio signal. We added therefore some fake y-ray data to repro¬ 
duce the radio signal at MJD ~ 54700. 

We also needed to change the amplitude of the response from 
flare to flare to match correctly the two peaks in the radio band 
at MJD ~ 55000 and ~ 55400. To obtain a qualitatively good 
match between our synthetic light curve and the radio data we 
used four different response normalizations. There is no phys¬ 
ical constraint for these responses to be the same for different 
radio outbursts: the physical details of the shocks are inside the 
response function, and different shocks may undergo different 
physical conditions. Table shows the scaling factors used to 
modify the original response and the time intervals. The first 
response is used to build the synthetic lightcurve before the be¬ 
ginning of LAT data (relying on the fake y-ray data introduced). 
The scaling factor is arbitrary chosen to be relative to the ampli¬ 
tude of the response in the first time interval. 
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Table 1. Time start and relative amplitude scaling of the several re¬ 
sponse profiles used to construct the synthetic lightcurve (red line in 
Fig.g. See text for the details. 


Time interval in MJD 

Scaling factor 

53971 -54710 

1 

54710-54890 

2.20 

54890 - 55340 

1.04 

after 55340 

1.83 


With these prescription we are able to qualitatively repro¬ 
duce the radio light curve after MJD ~ 54700. Figure shows 
the radio lightcurve (black line) and the curve derived from the 
LAT lightcurve with the described technique (red line). The LAT 
lightcurve is also drawn for comparison. Since the maximum of 
the response is at -1-160 days (Fig. [^, the average lag between 
y-ray and radio flares is therefore of the order of 140 + 20 days. 


4. Spectral Analysis 

To differentiate and further investigate the nature of the X-rays 
and the hard X-ray to y-ray components suggested by the cross¬ 
correlation analysis (Sect. we built the SED of 3C 273 in 
two different spectral states. We defined the PCA flaring state 
as the time intervals when the source flux in the RXTE-PCA 
lightcurve exceeds 18 photons s ' cm“^ and remain larger than 
this value for at least 9 days. Similarly we define the LAT flaring 
state as the time intervals when the source flux of the Eermi-LAT 
lightcurve goes beyond 1.2 ■ 10 ® photons s * cm“^ and remain 
larger than this value for at least 3 days. In case of two consec¬ 
utive flares we expanded the selected interval to the full width 
of the flaring period. The epochs defining the flaring states are 
shadowed in Eig. [T] 

Eor both spectral states, we built a multiwavelength spectrum 
from the X-rays to the y rays, using data from RXTE, INTE¬ 
GRAL, and Eermi selected in the time intervals identified as the 
PCA and LAT flaring states. 

The signal-to-noise ratios are enough to build multiwave¬ 
length spectra for each of the datasets, but not to follow the tem¬ 
poral evolution of the spectral parameters during the individual 
flares. Moreover, INTEGRAL data are not available for the total 
duration of the flaring times identified. 

In the following discussion we will also compare our datasets 
with the SED obtained from the multiwavelength campaign on 
3C 273 performed in 1991 (|Lichti et al.|1995]l and built with data 
from GINGA (2 - 30 keV) and data from OSSE, COMPTEL and 
EGRET on board CGRO (20 keV - 10 GeV). X- r ay and y-ray 
data showed a break at ~ 1 MeV. [Lichti et al.| ( |1995| l tested 
several models on their whole SED which includes also radio, 
infrared, optical and ultraviolet data. 

We first tested a simple cutoff powerlaw model N{E) = 
NcpE~^ !Ec), which is not able to fit the data: x]-ed 

is ~ 48 for the LAT flaring state, ~ 37 for the PCA flaring state 
and ~ 4 for the CGRO dataset. In the LAT and PCA flaring states 
the flux and the curvature of the LAT data is clearly not consis¬ 
tent with an exponential cutoff, and also in the CGRO dataset 
EGRET data are clearly above the exponential cutoff. 

We then used several models to fit the datasets: a broken 
poweriaw, a cutojf poweriaw and a iog-paraboia. None of these 
models has been used in|Lichti et al.|(fT995]l, where the high en- 
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broken powerlaw model 




Fig. 5. 3C 273 multiwavelength spectra at dilferent epochs fitted with 
the broken powerlaw model. The upper panel shows the SED and the 
best fit models (continuous lines), the lower panel shows the residuals 
from the best fit models. 


ergy data were fitted with an empirical function to model the 
break at ~ 1 MeV. 

The log-parabola model (see sect. |4.2| l is able to reproduce 
the typical inverse-Compton bumps measured by Fermi in 
blazar spectra, which show an intrinsic curvature ( |Nolan et al.| 
2012| l and can be interpreted with stochastic acceleration of 
the electrons ( [Massaro et al.||20()4l [Tramacere et al.|[2011| and 
reference therein). 

Figures |^|6 and show the multiwavelength spectra for the 
LAT and PCA flaring states, plus the CGRO dataset from |Lichd| 
jet al.j ( |1995[ l. Tables and list the best fit parameters, re¬ 
duced ^^”and degree of freedom for each dataset and each model 
tested. Errors are at 2.7 cr and the parameters without error are 
fixed. 

In general the CGRO dataset is closer to the PCA flaring state 
whatever model is used to fit the data. The main difference is the 
overall luminosity, which is systematically lower than during the 
LAT and PCA flaring states. We have also built the spectrum of 
the quiescent state, i.e. when the source is not flaring either in 
RXTE-PCA or in Fermi-LAT. It has a similar shape and is 2 
times brighter than observed by CGRO. 


4.1. Broken powerlaw model 

Fig. shows the three datasets with the best fit for a broken 
powerlaw model: 


A(£) = 


Nbp E-^' 
Nbp 


at £ < Ei 
at E > Eb 


The best fit parameters are reported in Table The fit by 
means of the broken powerlaw model does not provide a good 
description of the LAT and PCA flaring states spectra, but the 


Fig. 6. Same as Fig. [^for the log-parabola model. 


cutoff powerlaw + log-parabola model 



Energy (keV) 

Fig. 7. Same as Fig. |^for the cutoff powerlaw -l- log-parabola model. 
The dashed lines shows the two components of the model for the LAT 
and PCA flaring states. 


variation of the spectral indices between the two datasets gives 
some hints on the physics. This difference leads to a variation of 
the spectral break AF = F 2 - Fi from ~ 0.7 + 0.1 for the PCA 
flaring state to ~ 1.0+0.05 for the LAT flaring state, which could 
be related to a variation of the average slopes of the electron 
distribution (assuming the electron distribution modelled with a 
broken powerlaw too) from the LAT to the PCA flaring state. 
The break energy Eb is larger for the LAT flaring state than for 
the PCA flaring state. 
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Both datasets show a bump at ~ 20 keV (in RXTE and JEMX 
data), clearly visible in the lower panel of Eig. to understand 
if this structure is real, we built the spectrum of the Crab nebula 
to compare the deviations from the best fit model for the two 
sources. As the deviations are of the same order, we conclude 
that this bump is likely due to a calibration issue. 


4.2. Log-parabola model 

The log-parabola model is given by 


N(E) = Nip 



( 1 ) 


where Nip is the photon flux at £ = Eq, the a index denotes the 
slope at £ = £o and the p index models the curvature of the 
log-parabola. 

Eigure shows the three datasets with the best fit obtained 
with the log-parabola model. The best fit parameters are reported 
in Table|^ Eor the EAT flaring state the log-parabola alone is not 
able to fit the data adequately > 2). The residuals indicate 
that the intrinsic curvature of the EAT data does not correspond 
to that of the complete SED. 

The log-parabola model provides a good fit for the PCA flar¬ 
ing state, but it has to be pointed out that the EAT data in this case 
have worse statistics than for the EAT flaring state. Consisting 
of 3 points only, they do not show a clear intrinsic curvature. In 
this case the yS parameter can be determined only from the multi¬ 
wavelength fit and it is impossible to state if an extra component 
is needed or not. A similar statement can be done for the CGRO 
dataset, where EGRET data do not have enough statistics to ac¬ 
curately constrain the shape of the spectrum at y-ray energies. 


4.3. Cutoff powerlaw -h log-parabola model 

Erom the cross-correlation analysis (Sect. HD we inferred that 
X-rays and y rays are dominated by different components. In¬ 
terpreting the EAT data emitted by inverse-Compton processes 
in the jet, we tried to add an extra component in the form of a 
cutoff powerlaw at lower energies to model the X-ray emission 
in a similar way as assumed in Seyfert galaxies. 

Eigure |7] shows the three datasets with the best fit for the 
cutoff powerlaw -H log-parabola model. The best fit parameters 
are reported in Table Adding the cutoff powerlaw component 
improves the quality of the fit for the EAT flaring state, support¬ 
ing the hypothesis of an extra X-ray component. The parameters 
of the cutoff powerlaw component have been fixed to the typical 
values for Seyfert Galaxies ( |Mushotzky|1984[[Turner & Pounds] 
[1989^ and we used the same values of spectral index and cutoff 
energy for the three datasets because fitting the datasets with six 
free parameters leads to a degeneracy, especially for the cutoff 
energy This model provides a good fit also for the PCA 
flaring state and for the CGRO dataset. 

In order to constrain the best value for the cutoff powerlaw 
component we also tried to fit the X-ray data only (PCA, 
JEM-X, ISGRI and SPI) for the EAT flaring state with the cutoff 
powerlaw -H log-parabola model. In this case, the parameters of 
the log-parabola component are fixed to the values that provide 
the best fit of EAT-Eermi data alone. The idea is to model the 
contribution of the blazar-like component fitting the y-ray data 
alone and extrapolating its contribution at X-rays. Eitting the 
two components independently allows us to derive the best 
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Confidence Contours; Chi squared 



Phoindex 


Fig. 8. Dependency of they^ on the photon index T and cutoff energy 
El- of the cutoff powerlaw component for the PCA flaring state in the 
cutoff powerlaw + log-parabola model. The cross represents the best 
fit model. The lines from the internal to the external show confidence 
contour at cr = 1, 2, 2.706 (90%), 5. 


fit parameters for the Seyfert-like component, modelled as a 
cutoff powerlaw (Table |^. It is not possible to follow the same 
procedure for the PCA flaring state due to lack of statistics of 
the EAT data. 

We can also fit the whole multiwavelength spectrum using 
the spectral parameters of Table for the EAT flaring state in¬ 
stead of the fixed parameters previously chosen (Table|4 1 . In this 
way we obtain a new set of best-fit parameters and reduced 
for the cutoff powerlaw H- log-parabola model (Table |^. As ex¬ 
pected the is better for the EAT flaring state, but it is worse 
for the PCA flaring state, suggesting that the cutoff powerlaw 
component may not have the same parameters in both datasets. 
The last line of Table [^reports the best fit obtained leaving the 
cutoff powerlaw photon index T free to vary. 

In the PCA flaring state the variation of the spectral index T 
plays an important role: Pig. shows the dependency of they^ 
on the photon index E and the cutoff energy Ec. A reasonable 
good fit in the PCA flaring state can be obtained only if E > 
1.6 . At Ec = 360 keV the minimum is located at E = 1.66 
with xlcd - i -02, but the fit remains good for higher values of 
the photon index. The dependency on the cutoff energy is very 
weak: y^ is almost constant for a large range of Ec values. The 
data basically constrain only a lower limit to Ec. 


5. Discussion 

The analysis performed on the lightcurves and energy spectra 
of 3C 273 suggests different origin for the X-ray emission with 
respect to the y-ray emission. The y-ray emission at £ > 100 
MeV is compatible with the relativistically beamed SSC com¬ 
ing from the jet. This blazar-like origin is suggested by the 
cross-correlation between radio and y rays. The X-ray compo¬ 
nent is compatible with thermal photon emitted by the accretion 
disk and scattered in the X-ray domain by inverse-Compton pro¬ 
cesses, possibly on a hot corona close to the central black hole, 
as inferred for Seyfert galaxies. 
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Table 2. Best fit parameters, reduced and degree of freedom for the broken powerlaw model. The first column identifies the dataset. See text 
for details. 


Spectrum Ti Eb(MeV) r 2 Nbp xled 

LATf.s. 1.55 + 0.02 (7.8!|j) 2.61 +0.04 (2.641°;“)10-2 1.376 40 

PCAf.s. 1.65 + 0.02 (1.6!^;“) 2.41+0;“ (3.80 + 0.02)10-2 1.987 34 

CGRO 1.62 + 0.01 (2.2+2-9) 2.43+o}2 (0.91+0.03)10-2 0.718 16 


Table 3. Best fit parameters, reduced and degree of freedom for the log-parabola model. The first column identifies the dataset. See text for 
details. 


Spectrum Eq (keV) a p Nip xled 

LATf.s. 1 1.24 + 0.03 0.107 + 0.004 (1.58+°[2)10-2 2.311 41 

PCAf.s. 1 1.41 +0.02 (9.74+°;25)10-2 (2.901°;20)10-2 1.060 35 

CGRO 1 1.38 + 0.03 (9.8l!°;“)10-2 (0.70 + 0.03)10-2 1.290 17 


Table 4. Best fit parameters, reduced x^ and degree of freedom for the cutoff powerlaw + log-parabola model. The first column identifies the 
dataset. See text for details. 


Spectrum 

F 

Ec (keV) 

Ncp 

Eo (keV) 

a 

p 

Nip 

xld 

d.o.f. 

LAT f.s. 

1.7 

150 

0.028 

1 

0 75+0-06 

0.16 + 0.01 

(2.0 + 0.4)10-2 

0.971 

41 

PCA f.s. 

1.7 

150 

0.028 

1 

^■^^- 0.07 

0.13 + 0.01 

(6.3 + 1.0)10-2 

1.266 

35 

CGRO 

1.7 

150 

0.009 

1 

0 61+0-20 

0.19 + 0.03 

(0.34!°;28)10-2 

0.098 

16 


At a first attempt, the Seyfert-like component has been mod¬ 
elled with the identical cutoff-poweiTaw component in the LAT 
and PCA flaring states, using the same parameters T, Ec and Ncp 
(Table 1^. This means that we are able to explain the X-ray spec¬ 
tral variability with the variability of the blazar-like component 
only. However it is unlikely that the Seyfert-like component is 
not varying. In fact we can find slightly better fits using a dif¬ 
ferent spectral index T in the cutoff component (Table |^. Fur¬ 
thermore, the CGRO dataset clearly shows that at the time of the 
multiwavelength campaign described in |Lichti et al.| ( |1995[ ) the 
source was globally in a much lower state. 


We note that [Beaklini & Abraham| ( 2014| l reported a cross¬ 
correlation between 3C 273 Fermi-LAT data and radio data at 
7 mm wavelength (43 GHz) claiming a delay between 120 and 
170 days fully consistent with our estimate of ~ 140 d (Sect. 
|3.2| i, and with similar conclusions on the link between radio and 
y-rays. They explain the long-term variability of 3C 273 by vari¬ 
ation of the Doppler factor in a 16-year jet-precessing model. 
The LAT flaring state lasts roughly ~ 100 days, and according 
to the interpretation of Beaklini & Abraham ( |2014) l it is at the 
time when the Doppler factor is highest. The PCA flaring state 
spans a longer and not continuous time period: The simultane¬ 
ous LAT data may thus be a composition of photons emitted with 
a noticeably different Doppler factor. Our data are thus broadly 
consistent with the precessing jet interpretation of these authors, 
although we cannot really test it with the chosen two different 
states. 


5.1. The convolution and the 1C response 


Radio emission in radio loud AGN as 3C 273 is known to orig¬ 
inate through synchrotron mechanism due to shocks in a rela¬ 
tivistic jet (BegeImanjetar|T9M Marscher & Gear|1985| l. The 
response function used in the convolution (Sect. |3.2| l hides the 
processes that happen in the propagation of the shock, like the 
cooling of the electrons and the IC mechanism. In the inner part 
of the jet the electrons are energetic enough to produce a flare 
in y rays through IC processes. Further out the flare is produced 
at longer wavelength because of the electron cooling and the ex¬ 
pansion of the emitting region. 


The relative success of our convolution attempt suggests that 
the radio outbursts are closely related to the y-ray activity. The 
effect of the convolution is to smear out the complex variability 
structure of the Fermi lightcurve, resulting in only two appar¬ 
ently distinct radio outbursts. Our convolution thus implies that 
the radio outbursts would actually be a blend of emission from 
electron accelerated in several shorter shocks. 


It has to be pointed out however that longer flares in radio 
band may be the result of the superimposition of shorter flares: 
[Valtaoja et al.| ( |1999 1 decomposed radio lightcurves (22 and 37 
GHz band) of an AGN sample into exponential flares. |Hovatta ^ 
al. ( 2009|l used this technique to derive the variability of Doppler 


boosting and Lorentz factor of the jets from the time scales of the 
flares. These works show that it is possible to decompose a long 
radio flare into several sequential components. VLBI data would 
allow us to identify the different radio components and then it is 


possible to look for direct correlation with y-ray data. Jorstad et 
[aTjpOOlj l used this method to look for correlation between radio 
and y-ray properties in a sample of blazars. 


Article number, page 7 of 1 11 1 

























































Table 5. Best fit parameters, reduced^- and degree of freedom for the cutoff powerlaw + log-parabola model at LAT flaring time. Here the fit 
is done using PCA, JEM-X, ISGRI and SPI data. Parameters of the log-parabola are obtained fitting LAT data alone and used here as frozen 
parameters. See text for details. 


F 

EffkeV) 

Ncp 

Eq (MeV) 

a p 

Nip 


d.o.f. 

1.55 + 0.04 

(3.7;3-0)102 

(2-52!°i^)10-2 

100 

2.33 0.19 

2.63 10-“ 

0.934 

35 


Table 6. Best fit parameters, reduced^" and degree of freedom for the cutoff powerlaw + log-parabola model, using for the cutoff powerlaw the 
best fit parameters of Table Third line reports the best fit for the PCA flaring state if T is not frozen. The first column identifies the dataset. See 
text for details. 


Spectrum 

F 

Ec (keV) 

Ncp 

Eq (keV) 

a 

P 

Nip 


d.o.f. 

LATf.s. 

1.55 

370 

0.0252 

1 

0 

^**^^-0.08 

0.19 + 0.01 

(L6!°;9)10-4 

0.859 

41 

PCA f s. 

1.55 

370 

0.0252 

1 

1 7144-0.20 
^*'^^-0.08 

0.08!°:°^ 

(3.5!2-9)10-3 

2.057 

35 

PCA f s. 

1 

^•^^*0.05 

370 

0.0252 

1 

1 974-0.08 
^ ' -0.06 

0.11 + 0.01 

(8.0;|;^)10-3 

1.021 

34 


When building the synthetic radio lightcurve we had to vary 
the amplitude of the IC response from flare to flare (Table [^. 
There are several physical parameters which can affect the 1C 
process and may vary from shock to shock producing the ob¬ 
served variation in the amplitude of the IC response, but prob¬ 
ably the most important are the size of the shocked region and 
the injected electron distribution. The bigger is the size of the 
shock or the more energetic are the injected electrons, the larger 
is the number of scattered photons, and hence the more promi¬ 
nent is the flare in y rays. Different knots of the jet can be ejected 
at slightly different angles with respect to the observer’s line of 
sight or with a slightly different Lorentz bulk factor, resulting 
in a different Doppler factor (Jorstad et al. [2001 | l. The emis¬ 
sion of the knots with higher Doppler factor is hence boosted 
and shifted to higher frequencies, increasing consequently the 
total y-ray emission. Both these considerations can explain the 
observed variation of the IC response amplitude. 


5.2. Interpretation of the log-parabola model 

The energy spectrum of the LAT and PCA flaring states can be 
fitted with different parameter values. At first we note that the 
spectral index T of the cutoff powerlaw and the spectral index a 
of the log-parabola are both harder in the LAT flaring state. The 
best fit of the LAT flaring state is obtained with T = 1.55 and 
a = 0.39 while the best fit of the PCA flaring state is with T = 
1.66 and a - 1.27 (first and third lines of Table|^ respectively). 

In the PCA flaring state, the log-parabola model, because 
of the low statistics of the Fermi-LAT data, is sufficient to fit 
the dataset. Using the cutoff powerlaw H- log-parabola model a 
degeneracy is added and the parameters can vary in a quite wide 
range. 

The best fit yields a slightly steeper photon index F of the 
cutoff component together with a lower photon index a of the 
log-parabola and a slightly increased curvature ff. This means 
that with a softer cutoff powerlaw component the log-parabola 
index becomes harder in order to adjust the global slope of the 
spectrum (Table [^. However exploring the range of the possible 
spectral parameter in the PCA flaring state the parameter is al¬ 
ways smaller and the a parameter is always larger with respect 
to the LAT flaring state, meaning that the spectral variation be¬ 
tween the two states is real. It suggests a harder spectrum with 
stronger curvature in the LAT flaring state. 
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In order to find a possible physical interpretation of this dif¬ 
ference we note that it is possible to show that a log-parabolic 
distribution of particles n(y) generates a log-parabolic photon 
spectrum using the d-assumption for the emission of a single 
particle. In this case the relation between the spectral indices 
is a = (s - l)/2 and (3 = r/4, where s and r are respectively 
the spectral index and the curvature in the log-parabolic particle 
distribution n(y). 

Interpreting the log-parabola spectrum as the result of 
stochastic acceleration of the emitting electrons in a simple 
statistical approach, Tramacere et al. (20111 showed that s cx 
(it can be obtained from their Eq. 9) and 
where n is the mean number of acceleration 


• log(^ (o-g/e) ^ 
x-l 


(n {(Telef) 


steps, and e and cr^ are respectively the mean and the variance of 
the energy gain in each acceleration step. The spectral parame¬ 
ters are therefore related to the acceleration mechanism through 
the coefficient of variation crg/e. The larger is iXg/e, the softer 
and less curved is the spectrum. 

Our data show that the log-parabola component of the LAT 
flaring state is harder and more curved than in the PCA flaring 
state, and also that the peak energy of the SED is higher in the 
LAT flaring state. In the stochastic acceleration interpretation, 
this means that the acceleration mechanism during the y-ray 
flares is characterized by a lower coefficient of variation. This 
could result, for istance, because particles are uniformly boosted 
to higher energies through shocks. 


5 . 3 . Synchrotron Self Compton variability 


The SSC mechanism predicts constrains on the spectral variabil¬ 
ity which can be used to test the compatibility between the data 
and the model (i.e. see [Jones et aL]|1974[ [Bloom & Marscher] 
1996 and reference therein for a detailed discussion on SSC). 
Bloom & Marsch^ ( |1996| l explored the consequences of the 
variations of the physical parameters of the emitting source (i.e. 
the magnetic field, the particle density distribution, the viewing 
angle) on the resulting spectrum in the frame of the SSC model. 

We try to derive some informations on the physics of the 
X-rays and y-rays emitting region in 3C 273 based on the differ¬ 
ences between the LAT and PCA flaring states assuming that the 
blazar-like emission modelled as a log-parabola is due entirely to 
the SSC mechanism, without any External Compton contribute. 
In order to do this we introduce first Eqs. |^and[^which give the 
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value of the log-parabola peak energy Ep and the SED flux at the 
peak Ep (Massaro et al.|2004 1 : 


Ep = £o 10^2'“^''^^ 


( 2 ) 


log [EpFE^Ep)) = log Nip + . (3) 

The peak energy can be used to estimate the variation of the 
Lorentz factor y of the emitting electrons. In the SSC the ra¬ 
diation field is produced by synchrotron emission, so that the 
energy released by a single electron is peaking at a frequency of 
9 eB 

Vc oc y -. This radiation field then interacts with the electrons 

irieC 

themself, and photons are scattered to higher energies through IC 
scattering, increasing their frequency by another factor y^. 

The ratio between the energy peaks of LAT and PCA flaring 
states can provide information about the variation of y in the two 
spectral states. This ratio can be written as 


Ej jjEj (yL^BE 
E^p ylEl \r/>/ Bp’ 


(4) 


where the subscripts L and P stand for LAT and PCA flaring 
states, respectively, and the superscripts C and S stand for IC 
and synchrotron. Taking the energy peak of the log-parabola 
as E^ and hence using Eq. and the best fit parameters of the 
spectral states in Table we get yilyp - 1-6 if we assume 
Ep = Ep or yilyp - 1.2 if we assume Bi - Bp. Particles would 
thus be accelerated to higher energies during the LAT flaring 
state. 


To further investigate the physical parameters of the emitting 
region we can use the flux at the peak. We first rewrite the peak 

hv 

energy of the IC spectra in unit of e = -r as 


\ + z 


4 

r fip 


(5) 


where ep - , and i5o is the Doppler factor. The y factor here 

is that of the electrons emitting at the peak energy. 

This photon energy is used in Eq. wh ich gives the ob¬ 
served SSC vFy spectrum (|Dermer & Menon||20091 Ch. 7, eq. 
7.97) 


EFiie) 


cctIR'UbV 

llndp 



f 


'min(e',l/6') 


de'e' 




( 6 ) 


Equation]^ is obtained in the 5-function approximation and 
limiting the IC interation to the Thompson regime (upper limit 
of the integral). Prime quantities are measured in the rest frame 
of the source. B' is the size of the source (V' ~ R'^ being the vol- 
B^ 

ume), Up - — is the energy density of the magnetic field. The 
Stt 

electron distribution rieiy) appears twice, the former to account 



Fig. 9. Trend of variation of the SSC flux in the Ep - EpFpiEp) plane. 
The black points correspond to the LAT and PCA flaring states. The red 
line connects the two points with different Lorentz factor y, the blue line 
connects the two points with different magnetic field B. The parameters 
are listed in Table The dotted green line shows the trend associated 
to do variation. See text for details. 


for the IC scattering, the latter for the synchrotron emission, e' 
is the energy of the synchrotron photons which are scattered. 

Equations and relate the peak energy Ep and the vFy 
spectrum to the physical quantities of the source. Varying the 
physical parameters we can explore how the peak energy and the 
vFy flux calculated at the peak energy vary thereby and compare 
these trends with the data. 

According to Eq. [^the observed peak energy depends on the 
Lorentz factor y of the electrons at the peak and the magnetic 
field B (we do not consider here any variation of the Doppler 
factor 6d of the emitting region). The value of the observed Ep 
in LAT and PCA flaring states is reproduced by fixing arbitrary 
one of the two parameters and setting properly the remaining 
one. Then we can compute the SSC vFy flux at the peak energy 
with Eq. 1^ provided that the remaining parameters are fixed. 

To solve Eq. [^we assume a log-parabolic electron distribu¬ 
tion, for which we can derive the parameters s and r from the 
spectral fit as mentioned in Sect. |5.2[ and thus solve the integral 
numerically. We used the parameters of the best fit from Table 
l^for the LAT and PCA flaring states to derive the parameters of 
the electron distribution. 

Eigure shows these points together with the observations 
of the LAT and PCA flaring states, Table |^list the exact energy 
Ep and associated SED flux, together with the corresponding 
values of the magnetic field B and Lorentz factor y. 


We choose to investigate variation in B and y because in 
our simple description, Ep depends on these two physical pa- 


rameters on ly. 6d and R' have been fixed to 6n = 9 (Jorstad 
et al. 12005 i and R' - 1.6 10^^ cm (Rani et al. 20131, the re¬ 


maining parameter «' has been calculated to reproduce the ob¬ 
served flux. The total electron density n' is the same in both 
states, i.e. the integral of the log-parabolic distribution n[iy) is 
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Table 7. List of the varying parameters y and B used to draw the points 
in Fig. 1^ The remaining parameters are fixed to the following values: 
6o = 9, R' = 1.6 10'^ cm, n' = 2.7 10^ electrons/cm^, z = 0.158, 
dt = 2.30 10^’ cm, estimated for the redshift of 3C 273 in a flat universe 
(O/v = 1 - SIm) with Hq-1\ km/s Mpc and = 0.27. 


B (Gauss) 

r 

Ep (keV) 

EpFE{Ep)[^;) 

1.4 

2.8 10^ 

00 

O 

O 

1 


1.4 

2.2 10^ 

2.9 10^ 

1.5+2-3 lO-'O 

2.3 

2.5 10^ 

00 

O 

o 

7.9+fg° 10-® 

0.8 

2.5 10^ 

2.9 10^ 

3.11^3° 10-" 


forced to be constant adjusting the normalization constant. The 
electron density value needed to reproduce the observed flux is 
n' = 2.7 10® electrons/cm^. This looks like a big value, but 
is low enough for an optically thin emitting region: the optical 
depth T for inverse-Compton scattering in Thompson regime can 
be estimated to be t ~ R'crjn'^ ~ 10“^. 

Figurej^shows that the observed variability between the LAT 
and PCA flaring states is better described by a variable Lorentz 
factor y with a constant magnetic field B, rather than a constant y 
and a variable B. The confidence intervals are obtained account¬ 
ing for the errors on the spectral parameters: the parameters are 
modified by their errors at Icr and then the flux is recalculated 
with Eq. This indicates that the variation of the Lorentz factor 
y of the electrons is likely responsible for the difference between 
the LAT and the PCA flaring states. We note that variations of 
R' or n' do not produce a variation in Ep. 

We investigated also if the variation could be explained by 
the variation of the Doppler factor bo of the emitting region. 
Different knots may have different Doppler factors because they 
may be emitted at a different angle or with different velocity, 
as assumed in the varying bo-factor model of |Beaklini & Abr~ 
|ham| ( |20T^ . To test this, we used the same electron distribution 
n[iy) for the LAT and PCA flaring states assuming that the only 
varying parameter between the two states is bo. In this case Ep 
increases linearly with bo, while EpFEiEp) increases as b^ (see 
Eq. 1^ the trend shown by the dotted green line in Lig. |^. This 
trend is not compatible with the data. 

In conclusion, a variation of the electron distribution and 
consequently a variation of the electron Lorentz factor y cor¬ 
responding to the peak of the electron energy distribution seems 
to be the most likely driver for the difference between the LAT 
and PCA flaring states. 

5.4. Location of the y-ray emitting region 

The spectrum observed by Lermi-LAT is not significantly ab¬ 
sorbed up to ~ 3 GeV. Photons below 3 GeV could produce 
pairs by interacting with target photons of energies above 0.3 
keV. The latter are probably emitted from the accretion disk and 
corona within 100 gravitational radii from the central black hole. 
The y-ray emitting region should be located outside of that re¬ 
gion, which is not very constraining. 

Since the sensitivity of Cherenkov telescopes is not sufficient 
to detect 3C 273 in the TeV domain, where the pair absorption 
signature should be present if the y-ray emitting region is inside 
the broad line region, the position of the y-ray source cannot 
be constrained with respect to the size of the broad line region 
( |Donea & Protheroe|2003| l. 


6. Summary and Conclusions 


The high energy spectral variability of 3C 273 has been stud¬ 
ied from ~ 1 keV up to ~ 10 GeV and compared with ra¬ 
dio data. We presented Lermi-LAT, RXTE-PCA and Radio (37 
GHz) lightcurves and multiwavelength spectra of RXTE-PCA, 
Lermi-LAT and JEMX, ISGRI, SPI on board INTEGRAL, built 
for two different spectral states, namely the LAT flaring state, 
i.e. when the source is flaring in Lermi-LAT, and the PCA flar¬ 
ing state, i.e. when the source is flaring in RXTE-PCA. 

Timing analysis (Sect. indicates the presence of two dif¬ 
ferent components in the considered energy band. Gamma ray 
data correlate with radio data, but not with X-ray data, confirm¬ 
ing that the emissions in y rays and radio have a common origin, 
while the origin of the X-ray emission is different. The convo¬ 
lution between the y-ray lightcurve and the SSC response func¬ 
tion reproduces the radio signal (Sect. |3.2| l. This indicates that 
the data are compatible with the synchrotron mechanism in the 
frame of a shock-in-jet scenario ( |Marscher & Gear|1985[[Ttirler| 
et al.|1999 20001. The necessity to shrink the response function 
obtained by low energy data suggests that what is seen as a sin¬ 
gle flare in the radio lightcurve may be in this case the blended 
emission from several shorter shocks, which have been detected 
separately in the y-ray lightcurve. 

Spectral fitting of the high energy data (Sect. shows that 
a single model emission as a cutoff powerlaw or a log-parabola 
cannot fit adequately the whole SED from ~ 1 keV up to 
~ 10 GeV. Knowing from timing analysis that X-rays and y 
rays have different origin, we used a two components spectral 
model: a cutoff powerlaw at low energies and a log-parabola at 
high energies, the former modelling the Seyfert-like emission, 
the latter modelling the jet emission. We found that for the 
LAT flaring state the powerlaw and log-parabola components 
are harder than for the PCA flaring state, and the LAT flaring 
state spectrum has a stronger curvature at the SED peak. The 
spectral variability between the two states can be understood if 
it is driven by a change in the electron distribution with a peak 
energy increasing by ~ 30% during the LAT flaring state. 


Unlike typical blazars, the high energy spectrum of 3C 273 
is peaked in the MeV energy range, below the Lermi-LAT lower 
energy threshold. Even during the flare, the peak of the emis¬ 
sion remains undetected. As there are no experiments sensitive 
enough in the MeV, it is difficult to study in detail the spectral 
properties of the y-ray emission of 3C 273. Nevertheless, the 
shape of the y-ray spectrum of 3C 273 is very similar to the 
shape of the GeV peaked blazar detected by Lermi, which are 
usually fitted with the log-parabola model, suggesting that the 
electron distribution looks universal in very different conditions. 
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